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ABSTRACT 


A new  formulation  for  radiation  from  open-ended  guiding  structures 
is  presented.  The  differences  between  this  formulation  and  conventional 
physical  optics  or  diffraction  ray  theory  is  discussed.  The  theory  has 
quite  general  validity  and  is  here  applied  in  detail  to  open-ended  trans- 
mission lines,  coaxial  cables,  parallel  plate  waveguides  and  circular  wave- 
guides. For  the  examples  discussed,  checking  with  existing  solutions  is 
provided.  Finally,  the  theory  is  extended  to  the  case  of  flared  open- 


ended  guides. 


■I 


CHAPTER  I 
THE  PROBLEM 

Radiation  from  apertures  is  certainly  important  from  the  application's 
viewpoint.  In  this  paper  radiation  from  open-ended  guiding  structures, 
as  waveguides,  coaxial-cables  and  horns,  is  considered.  The  classical  way  of 
analysis  is  based  on  the  so-called  apCAtutiC  method  [l]  which  uses  the  phy- 
sical optics  approach  coupled  with  Kirchhoff's  approximation.  More  recently, 
optical  or  qu&bt-OpttcCLZ  techniques  have  been  applied,  decomposing  the  wave 
inside  the  guide  into  a set  of  plane  waves  and  then  considering  the  scattering 
of  these  plane  waves  by  the  edges  of  the  guide.  The  basic  idea  is  to  model 
the  latter  problem  as  the  scattering  of  a plane  wave  by  a metal  half-plane. 
These  ideas  and  tools  have  been  applied  to  the  analysis  of  reflection  inside 
[2-7],  and  radiation  from  the  guide  [8-17].  A different  model  has  also  been 
suggested  for  the  analysis  of  the  scattering  process:  a half-plane  over  either 
an  electric  or  a magnetic  infinite  wall  [l8-l9].  It.  is  also  important,  to 
recognize  that  exact  solutions  to  the  radiation  problem  are  available  in  a 
limited  number  of  cases:  the  parallel-plate  and  cylindrical  waveguides  [20], 
the  latter  with  either  azymutally  symmetric  or  asymmetric  excitation.  The 
scattering  by  two  half-planes  has  also  been  considered  [21]. 

The  aperture  method  has  well  known  limitations  and  shortcomings.  It  is 
reasonably  simple  only  when  the  scalar  approximation  is  adopted;  the  vector 
theory,  however,  yields  only  marginal  improvements  [ill. 

The  quasi-optieal  approach  seems  to  give  good  results.  At  transitional 
regions,  e.g.,  at  lit-shadow  boundaries,  the  expression  of  the  diffraction 
coefficient  becomes  rather  involved  [22,  23]-  In  the  caustic  regions,  either 
fictitious  equivalent,  edge  currents  [fh]  or  a spectral  theory  [88]  should  be 


adopted;  the  latter  is  also  useful 


if  not  necessary,  for  non-planar  exci- 


tation. In  the  case  of  multiple  reflections,  as  always  happens  for  apertures, 
the  analysis  although  straightforward  in  principle,  is  extremely  cumbersome 
in  practice. 

There  is,  however , at  least  another  possibility  for  computing  radiation 
from  open-ended  guiding  structures.  Probably  the  best  way  to  introduce  the 
basic  philosophy  of  the  new  method  is  to  consider  a wire  antenna  with  sinu- 
soidal current  distibution  I ^ sinS^z  (see  Fig.  l)  and  a time  dependence 
exp(iu)t).  The  far-field  phasors  can  be  easily  computed,  the  result  being: 

I sin0  ~ icosO  cosGsinS  1 + cosG  cosS  1 

E = C — 0 exp(-iBr)  — + exp(-iBcosO)  — — — (l.l) 

Hnr  cos2©  - cos2©  cos2©  - cos20 

a a 

h = 4>  , C ',  6 = u>{evT  (i.:: 


)« 


where  3 = BeosO  . 

a a 

Eq.  (l.l)  suggests  this  model  for  the  radiation  process:  the  far-field 
is  given  as  the  superposition  of  two  spherical  waves  from  the  end-points 
z = 0 ; - 1.  of  the  antenna.  For  a wire  infinitely  long  in  the  -z  direction 

the  second  term  in  the  bracket  of  (l.l)  can  be  neglected,  and  we  have: 


<E0  = ^ 2 7T 


sinOcosO 

a 

cosO  + COS0 


exp(-i3r ) 
COS0  - COS0 


(1.3) 


lH0  = E0^ 


(l.U) 


where  1^  = - IQ/2i  is  the  amplitude  of  the  -chC-tdmt  cuAAent.  Note  that  the 
electromagnetic  field  (1.3-M  can  be  conveniently  derived  from  a vector  po- 
tential 


1 COS0 


exp(-i(3r ) 


-C 


A 

Z 


-io)  2ur  cos0  + cos0  cos0  - cos© 

a a 


wherein  the  connection  with  far-fields  is  given  by 


. _ .A  A 

E = -it)  r x A x r 


u 10)  A 

H = - — r x A 

L, 


(1.5) 


(1.6) 

(l.T) 


For  a two-wire  transmission  line — the  simplest  conceivable  guiding  struc- 
ture— ,TEM assumption  for  the  incident  wave  (0  = 0),  hypothesis  of  no  higher 

cl 

order  interaction  between  the  two  wires,  and  superposition  give  (see  Fig.  2): 


exp(-igr ) 


2frr 


exp( i3qsin0cos<|> ) - exp(-iBasin©cos<{>) 
I - sin2Qsin2i}> 


A 

y 


Fig.  ?.  Radiation  from  a two-wire  transmission  line. 

Note  that  the  field  is  everywhere  well-behaved.  In  particular  the  (normal 
radiation  diagram  is  given  by: 

sin(  8acos(J> ) 

E-plane  (0  = 90°)  : 

Sacos<J> 

H-plane  (0=0)  : sin(BasinO) 

Thz  hAdjjxtXon  pnoczM  h>  modelled  cu>  two  iphztucal  wavzi  fi/iom  the  two  znd- 
potntb  oi  thz  tSLtzn&nuA&'ion  Linz.  This  model  has  a simple  physical  counter 
part  in  time-domain  when  a pulsed  voltage  is  propagating  along  the  transmi 
sion  line.  The  local  charges  induced  on  the  two  curves  (apparently)  move 


L 


6 


F ' 


I 


with  constant  velocity  and  radiation  takes  place  only  at  the  very  end  of 
idle  line  where  charges  are  reflected  with  a (di  scent  inui'iir)  change  of  velo- 
city  direction.  This  idea  has  been  exploited  in  [del.  Note  further,  for 
future  reference,  that  the  fields  individually  radiated  by  each  wire  di- 
verge in  the  direction  0 = 0,  i.e.  , the  direction  of  the  incident  wave. 
However,  due  to  the  dispersion  equation  of  the  line,  cos0  = 7,  the  two 

3, 

fields  compensate  and  produce  a finite  result  along  0=0. 

For  an  open-ended  guiding  structure  we  will  have  surface  currents  inci- 
dent on  the  truncated  rim  of  the  guide.  Accordingly,  the  Hadiatlcn  pHoceS-5 
can  be  modelled  at,  i>upcHpo!>ition  oft  elemetaHtj  iphenlcal  idoveA  coining  fifiom 
each,  point  ofa  the  Him.  Obviously,  a satisfactory  theory  should  take  into 
account  interaction  between  incident  surface  currents.  We  will  derive  ra- 
diation coefficients  associated  with  incident  surface  currents  from  the  sim- 
plest possible  canonical  problem  for  which  an  exact  solution  exists:  the 
parallel-plate  open-ended  waveguide . 

It  is  important  to  recognize  the  profound  difference  between  this  ap- 
proach and  the  quasi-optical  one.  The  latter  emphasizing  the  geometrical 
relationship  between  field-point  and  rim-point,  so  that  the  radiated  field 
is  due  to  a finite  number  of  contributions  along  the  rim  (except  in  caustic 
regions  where,  in  any  case,  the  quasi-optical  approach  fails  to  give  a satis- 
factory answer).  On  the  contrary,  this  new  approach  is  still  based  on  a 
field  theory  rather  than  a ray  theory  and  the  radiation  is  computed  as  an 
integral  along  the  truncated  rim.  Accordingly,  as  long  as  the  radiation 
coefficients  are  well-behaved,  no  problem  can  be  anticipated  for  the  compu- 
ted field  everywhere.  It  is  here  noted  that,  this  method  has  nothing  to  do 
with  the  equivalent  edge  current  procedure  [?U],  which  is  still  based  on  a 
ray-theory  model  and  makes  reference  to  fictitious  electric  and  magnetic  cur- 


' 


rents . 


CHAPTER  2 


RELEVANT  FORMULAS  AND  EXPANSIONS 


For  the  reader's  convenience,  most  of  the  mathematical  formulas  used 
in  the  next  sections  are  collected  below. 

Weinstein' 6 cU^racJMjn  function,  exact  [27]. 

Definition : 


U (a,  p)  = 


2TTi  J 


In [ 1 + exp(-ip  - pu2 ) ] 


(;  - 


1U 


( / - if  ) 


1/2 


u( / - if  ) 


1/2 


1 - i 


du 


Properties : 

U-  (-a,  p)  = -U-(a,  p) 


U— (+_  0,  p)  = + f ln[  1 + exp(-ip) 


Expansion  for  large  a and  p: 


U— (a,  p)  = 


|2up  a 


E 


— exp(  — iTr/U  ) X (+  )n  exp(-inp) 


n n 


3/1 


♦ 0 lo-’o-3'1 


Expansion  for  small  p: 


U+(o,  p)  = ]r  ln[f  + exp(-ip) ] + [ in  ^ + 1 - y - + 0 [po] 

d 2n  p 

U-(o,  p)  = ln[p(;  + a)]  + jp-  - [ In  ~ + 1 - Y - —■  ] + 0 [po] 


Weiru>  tetri' * dt^ractton  function,  approximated  [20,  27]. 


For  large  p the  main  contribution  to  the  integral  (2.1)  arises  from 


the  vicinity  of  the  saddle  point  u = 0.  Thus,  neglecting  all  terms  of  order 
u2  [27]: 


Ui  (a,  p)  = / ln[j  + exp  (-ip  - pu2 ) ] 

J U - — | — o 

— oo 


The  difference  between  U—  and  U—  is  of  0[p  7]  uniformly  over  the  range 
-J  < 0 _<  J. Letting  p = 2nq  and  s (and  using  the  same  symbol  U— ) : 


U^s,  q)  = 


/ 1""1 


exp(-27Tiq  -|  )]  (T-i')V 


which  is  the  expression  given  in  [20].  Tables  of  (the  complex  conjugate  of) 
U are  referred  to  in  [20,  28]. 

Properties : 

Same  as  (2.2-3) 

U— (s,  q + n)  = U— (s,  q)  I 

U+(s,  q)  = U (s,  q+~)  ( 

Expansion  for  large  s: 


U-(s,  q)  = 


exp(-iTT/U)  ^ (2m  - 1)!!  * exp(-2ninq) 


Ll^m  - i J ! ! exp  (,  -2 

r / . ,,  £ m m + 

m Lsexp(-nr/4r  n n n 


m + 3/2 


exp(  -iir/1*) 


\2tt  s 


Eexp(-2irinq) 

575 * 0 ls  1 


n n 


Expansion  for  small  s and  q: 


U— (s,  q)  = ^ In 2 + 0 [s] 


9 


II  (s,  q)  = - -(In?  - i^-)  + ]ri(<;  + rjlinq')  + 0 [ ;;  ] 

liiwkef  functions 


a (jc)  = arg[Hp^(x)]  = 


I 

8ttx  + 


Ti(x) = 4ar«[Hi2)^)i  + i=  ? + if  + sk~  ••• 


1 X 1 


a = 
m 


L rtt(2)/  n / r 2m  + 1 m2  - 1 /h 

- arglH  (x)  ] = - [x  - r — tt  + — — 


tt  m 


2x 


a - - V arg[Hi2)(*)]  *j=ju-  25Lii  , * ^-|x3/1* 


m tt  m 

Integral* 

, IT 


2i  sin(xcos4> )cos<})  d<}>  = 27TiJj(x) 


0 


(.M') 


(2.13) 

(2. lit) 

(2.15) 

(2.16) 


(2.17) 


(2.1S) 


2 / cos(xcos<j>)  d<t>  = 2ttJ  (x) 

./  o 

0 

m 2tt 

[exp(ixcosC)  + (-)  exp( -i xcos? ) ] exp(lmC)  d£  = exp( ixcos£ )exp( im£ )d£ 


2TTiraj  (x) 
m 


r tt  m 

/ [exp(ixcos^)  - (-)mexp(-ixcosf;)]  exp(  imf,  )cos£dC  = 


.277 

/ exp( ixcos£ )exp( imC )cos£d£ 


_ . 7 J . 7 ( x ) - J , (x) 
. m + / m + / m - 1 


-7T1 


(2.19) 


= -2TTim  + ? J'  (x) 
m 


(2.20) 


lu 


r 


[exp(ixcosO  - ( -)mexp(-ixcos  ]exp(  iiri 


2TT 


/J  I ( X ) + J f ( X ) 

exp(ixcos£)exp(im£)sin£d£;  = 2TTi  — — — ^ 


= 27Tira  — — 


mJ  (x ) 


J" cos(-^y2-  uy )dy  = f- 


2a 


2 .2 


COS' 


UX_  C (0-  _ H^)  + sin. 


. 2. .2 


a ir^ 


/* sin(— - 

J ^_2 


2a 


uy)dy  = - 


L 2nf/  ” 2it£j 

S (&L-  _ Hi)  + 


2 0 a tt£ 


(2.21) 

(2.22) 


cos- 


2 2 
tl  cl 


2tt£ 


sin  HJil  c (LiL_  _ H i) 

2 a iff;'  - . , 5 1 a iff/ 


2tt£' 


where  C(x)»  S(x)  are  the  Fresnel's  integrals. 
FsizAneX  in£e.gnal 6 and  auxituviy  function 


C(t)  = ± f^2  + f (*•)  sin— - g(t)  cos-^|-j 
S(t)  = ± [J-  - f(t)  cos^|—  - g(t)  sin— j 


(2.2U) 


(2.25) 


where  f(t)  and  g(t)  are  the  so  called  auxiLuxAij  function*,  and  plus  and  minus 
signs  apply  to  positive  and  negative  values  of  t respectively. 


(2.26) 


(2.27) 


(2.28) 


(2.29) 


(2.23) 


J 


lJ 
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CHAPTER  3 


WEINSTEIN'S  SOLUTION  FOR  THE  PARALLEL  PI  ATE  OPEN-ENDED  WAVEGUIDE 

The  problem  of  radiation  from  an  open-ended  parallel  plate  waveguide 
has  been  rigorously  solved  by  Weinstein  in  19^*8,  (results  are  referred  to 
in  [20],  Ch .7),  and  independently  by  Heins  [29],  essentially  by  applying 
Wiener-Hopf  techniques.  Solutions  to  this  problem  (and  also  to  the  case  of 
a circular  waveguide)  are  referred  to  in  a number  of  books  [30-33].  Here- 
after, Weinstein's  solution  is  summarized ^ + ^ with  reference  to  Fig.  3.  In 
addition,  only  ^OA-^dldit  and  -tt/2  _<  <f>  _<ir/2  {falont  Apace)  are  considered, 
where  the  latter  condition  can  easily  be  relaxed. 


Fig.  3.  Radiation  from  an  open-ended  parallel  plate  waveguide. 


(+)Note  that  in  the  original  Weinstein's  solution  [20l  (i)  the  assumed  time 
dependence  is  exp(-iist),  (ii)  the  z-axis  is  pointing  toward  the  waveguide, 
(iii)  the  unit  system  is  the  gnussian  one  and  (iv)  symbols  are  slightly 
different  from  those  here  adopted. 


Id 


For  an  incident,  TF,  wave  (no  x-variation)  the  field:'  can  be  derived, 
no 

via  (1.6-7),  from  a vector  potential  \ - Ax  where: 
for  n odd : 


CJ  exp[-i(3p  + it/1*)]  cos^-  cos( Basin*) 


A = 


2tt3p' 


"(J)  ~ 

cos  — — - cos 4>  - cos^ 


exp(U+  - U+) 


for  n even: 


iCJ  exp[-i(Bp  + tt/U  ) ] cos|-  sin (Basin*) 


A = 


0) 


\j  2tt3p' 


exp(U  - U ) 
n 


cos- — COS*  - cost}) 


In  (3.1-2)  the  angle  4*  is  related  to  the  waveguide  propagation  constant  Bn 


B = \/b2  - (mr/2a)2  = Bcost})  ; 

n V n 


J = J X in  the  in.cide.nt  bunAace  cuAAent  at  z = 0 related  to  the  incident 
— n n 

electric  field  at  z = 0,  E = E x,  by: 

-o  o 

E sin* 

, on  . 

J = , lower  rim, 

n » 

Jc 

(-)n+,E  sin* 

J = — , upper  rim, 

n K 

U— (s,  q)  is  the  Weinstein  diffraction  function  (2.7): 
lA(s,  q)  = U—  (>[2Bacos*,  ) , 


and 


(3.1) 

(3.2) 

by: 

(3.3) 

(3.*0 

(3.6) 

(3.6) 


IF-j-  (s,  q)  = U—  (^23acos*\  ^-) 


(3.7) 


Results  (3.1-2)  are  appA.0XAmcLte.cl  onty  because  the  approximate  ex- 
pression (2.7)  for  U—  has  been  used.  They  become  exact  if  the  exact  expres- 
sion (2.1)  is  adopted  with  0 = cos<J>,  p = 23a.  However,  ever  for  moderate 
value  of  Ba,  the  difference  between  the  two  formulations  is  negligible  save 
in  the  neighbors  of  |<J)|  = tt/2,  where  it  remains  small. 

For  an  incident  TM  wave,  the  fields  can  be  derived  from  a vector  po- 
tential  A = Az  where: 


for  n odd: 


SJ  exp[-i(8p  + ttA)]  cos  — cos(  Basing ) 


for  n even: 


exp(U+  - U+) 


cos  £ cosd>  - cosd> 

2 T n 


(3.8) 


i^J  exp[-i(Bp  + A/k)]  cos  sin(Basin<{>) 


■\f2i1ep' 


— I I 7 “p  » ' " ’ 

cos  jr  cos(|)  - cosd) 

2 n 


(3.9) 


where  U—  is  again  given  by  (3.6)  and  the  ineJAent  iuA^aee  cuAAent  = J^x 
at  z = 0 is  related  to  the  incident  magnetic  field  at  z = 0,  H = H^x,  by: 
Jn  = Hq  , upper  rim 


(3.10) 


J = (-)  H , lower  rim 
n 


(3.1D 


Solutions  (3.1,  2,  8,  9)  deserve  a comment.  For  Ba  large,  use  of  (2.10) 


shows  that: 


exp[U—  - U— ] = f + 0[  7/yBa] 


(3.12) 


\ccordingly,  the  left-hand  side  of  (3.12)  is  an  inteAaction  factoA  which 
approaches  unity  when  the  spacing  between  the  two  plates  is  large. 

Let  us  momentarily  neglect  this  factor.  It  is  immediately  recognized 


that  the  fields  are  proportional  to: 


f 


exp[-i(Bp  + tt  ' U ) ] cos~  exp(-iBp')  exp(-iSp") 


Y 2n 


PP 


<j> 

n . , 

•or.  — eosip  - cos  if 


(3.13) 


for  the  TE  case,  and  to 
exp[  -i  ( Bp  + tt  -'U ) ] 


. i|> 

cos  — sin  2 


exp(iBp’)  ± exp (-i Bp") 


( 3. lU ) 


2irBp  *■  *■  cos<)>  - cosi))^ 

for  the  TM  case.  In  (3.13-1*0,  the  plus  sign  applies  to  n odd  and  the  minus 
sign  to  n even.  The  total  field  is  then  obtained  as  the  superposition  of  two 
cylindrical  waves  coming  from  the  upper  and  lower  rim  of  the  guide  respectively. 
When  the  expression  for  the  incident  current  (3.^-5,  10-11)  is  taken  into 
account,  the  similarity  between  (3.13-1*0  and  the  metal  half-plane  scattering 
problem  [3*0  is  immediately  recognized.  Each  constituent  component  of  (3.13-1*0 
coincides  with  the  asymptotic  evaluation  of  the  field  scattered  by  a metal 

half-plane  when  the  incident  wave  is  plane,  at  an  angle  <(>  with  the  two  rims 

n 

respectively  and  with  polarization  E = x (TE  case)  or  II  = x'  (TM  case). 

This  is  obviously  consistent  with  the  splitting  of  the  incident  field  in  a 
set  of  two  plane  waves  at  angles  +4)^  with  respect  to  the  z-axis.  The  above 
similarity  has  been  recognized  from  the  beginning  [20]  and  is  the  basis  for 
all  quasi-optical  approaches  to  this  problem.  There  is,  however,  a profound 
difference  between  the  half-plane  problem  and  the  one  at  hand,  and  this  has  been 
so  far  apparently  overlooked.  Each  constituent  tcj un  in  ( 3. 13-1  *4 ) divcAgte 
at  4>  = <J>  (lit-shadow  and  lit-reflection  boundaries  respectively)  while  the 
totat  faictd  dote  not.  Consider,  for  instance,  (3.1).  When  4>  = 4>n  then 
Basing  = nn/2  and  is  the  normalized  eigenvalue  of  the  parallel  plate  guide. 
Accordingly, 

cos(  Basing)  , Bacoscb 

lim  = (i)n-J  B_ 

cosd)  - cos<j> 
n n 


sin<(> 


(3.13) 
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and  is  finite.  The  reason  for  that  is  very  basic  and  quite  general.  Consider 
the  sketch  of  Fig.  h where  ray  tracing  is  depicted  for  the  parallel  plate 
guide. 


Fig.  U.  Ray  tracing  inside  a parallel-plate  waveguide. 


Points  C and  E,  as  well  as  G and  F,  are  on  an  equiphase  surface;  accordingly, 
the  phase  delay  along  the  two  paths  GHIC  and  FE  should  differ  of  an  integer 
number  of  2n  (conb-Lbteney  KticUionit hip , the  optical  analog  of  field  dispersion 
equation).  Consider  now  the  phase  delay  along,  the  two  paths  CBP  and  EDAF.  The 
phase  introduced  by  the  scattering  process  at  A and  L does  not  depend  on  the 
angle,  so  that  it  is  equal  for  the  two  paths.  However,  the  usal  phase  delay 


1 1\ 


of  tt  at  A (due  to  the  reflection)  is  now  missed,  since  the  reflection  co- 
efficient is  substituted  by  a scattering  coefficient.  Accordingly,  the 
two  ray  contributions  at  P from  L and  A are  opposite  in  phase. 

The  following  conclusion  can  be  drawn.  For  the  half-plane  diffraction 
problem,  the  spectral  components  of  the  scattered  field  display  pole  singu- 
larity (in  the  wavenumber  space)  which  approach  the  saddle  point  position 
at  lit-shadow,  lit-reflection,  so  that  the  steepest  descent  method  evaluation 
(performed  in  the  usual  way)  becomes  invalid  (resort  should  be  made  to  the 
Pauli-Clemmow  modified  steepest  descent  method  of  integration  [35]).  On 
the  contrary,  these  pole  singularities  cancel  each  other  for  the  problem  at 
hand,  so  that  the  usual  asymptotic  form  for  the  scattered  field  is  valid 
through  all  space. 

The  interaction  factor  (3.12)  is  now  considered.  Use  of  (2.10),  valid 
for  26a  cosij)  large , shows  that: 


exp(U—  - U— ) = 1 + 

r n — 


exp[-i(2Ba  + tt  / h ) ] coscf  - COS(J) 


V ^ir6a 


COS(J)COS(f) 


(«)r 


exp( -i2nBa) 


(n  + 1) 


3/2 


+ 0 [ //da] 


(3.16) 


k 


When  (3.16)  is  substituted  into  (3.1-2,  8-9)  ail  extra  terns  appears  for  the 
vector  potential  A.  This  new  term  can  be  easily  constructed  using  the  model 
of  interaction  process  depicted  in  Fig.  5 (tingle  interaction) . The  ray 
along  BC  explains  the  term  n = 0 of  the  series  in  ( 3 . 1 6 ) ; and  the  remaining 
terms  of  the  series  are  due  to  the  nonuniform  illumination  of  the  edge.  For 
future  reference  it  is  important  to  recognized  that  2irq  = 28a  takes  care  of 
the  phase  delay  of  the  ray  path  between  the  two  edges. 

When  28a  is  small,  use  of  ( ? . 1»— 5 ) shows  that: 
exp  (U*  - U+)  1 (3.17) 


cos 


exp  (U~  - U ) = — 

n £ 


(3.18) 


cos 
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THE  ELEMENTARY  SCATTERING  COEFFICIENTS 


The  field  radiated  by  an  open-ended  parallel-plate  waveguide  has 
been  represented,  according  to  Weinstein's  theory,  as  the  superposition 
of  Cijlindtiacal  Wavdi  from  the  rims  of  the  guide.  The  main  idea  is  now  to 
represent  each  cylindrical  wave  as  superposition  of  e.Zo.n\tntaAy  i> pkeAtcal 
um>e^  radiated  at  each  point  along  the  rim.  This  goal  can  be  accomplished 
by  noting  that,  asymptotically  as  $P  ->  °°, 


-2uiF(0) 


exp[  -i  ( ftp  + tt/U  ) ] 


^2tt0p  ' 


+oo 

» j p(0 


exp(-ig  ^p2  + £ 2 ) 


2 + 52 


dC 


(4.1) 


Here  p is  identified  with  p'  (p")  in  (3.1-2,  8-9),  £ with  a curvilinear 
coordinate  along  the  upper  (lower)  rim  of  the  parallel-plate  waveguide, 

^~p2  + = r is  the  radial  distance  from  rim  point  to  field  point  and  F(0) 

can  be  immediately  obtained  by  inspection  from  (3.1-2,  8-9). 

The  next  step  is  the  determination  of  F(£),  which  should  coincide  with 
F(0)  times  an  arbitrary  factor  which  reduces  to  zero  as  ^*0.  This  factor, 
in  general,  will  not  play  an  important  role  when  applications  are  made  to 
ZjOAQZ  apertures.  As  a matter  of  fact,  superposition  of  elementary  radiation 
from  the  rim  of  the  aperture  will  result  in  an  integral,  and  in  the  asymptotic 
evaluation  of  such  integral  the  above  factor  is  unimportant.  Only  the 
(coscj)^  - cos<J>)  term  appearing  in  the  denominator  of  (3.1-2,  8-9)  should  be 
handled  with  care  to  avoid  a singular  behavior  of  the  elementary  scattering 


i 

i 


coefficients  which  are  going  to  be  introduced. 

As  a conclusion,  the  following  A Hi  at Z will  be  made. 


p 


_ 


20 


dJ„  exp(-i 0r ) 

d A = r. 

-iui  pur 


con  — exp(i0asin0sin*|O  +_  oxp(- i 0as  inOsin4> ) 


r — n~  ■>  < 

\|  / - S 1 n 'Or.  i tl ' i)'  - COiiO 


exp(U~  - U— ) 

s 


(4.2) 


dJ. 


dA  = C- 


-IU) 


exp(-i0r)  cos  -&■  exp(iBasin0sin4>)  exp( -i0asin0sin9 ) 


2nr 


cos  - 


fr-~  sin20sin2<J)'  - cosG 


exp  (U—  - U— ) 


(4.3) 


In  (4.2-3)  the  choice  of  plus  or  minus  sign  is  related  to  equal  or  opposite  sign 

for  the  incident  surface  current  at  Q',  Q"  respectively;  the  angle  0 is 

given  by  0 = 0cos0  , where  0 is  the  waveguide  propagation  constant.  Note 

g g g 

that  ?0asinOsinc(>  = jr_'  - r"  | and  that  the  elementary  vector  potentials  ( U . 2—3 ) 
are  everywhere  finite  provided  that  the  dispersion  relation 

exp(i0asinO  ) + exp( -i 0asinO  ) = 0 (4.4) 

g g 


is  valid. 

It  is  obvious  that  application  of  (4.2-3)  to  the  case  of  a parallel 
plate  waveguide  using  the  asymptotic  evaluation  (4.1)  of  the  superposition 
integral  gives  (3.1-2)  and  (3.8-9),  respectively. 


CHAPTER  S 


RADIATION  FROM  OPEN  ENDED  CIRCULAR  WAVEGU I DE. 


AZYMUT! j\1.IjY  SYMMETRIC  MODES 


Consider  now  an  open  ended  circular  waveguide  as  depicted  in  Fig.  7. 
The  guide  is  excited  "by  4>  — independent  modes  (fl.zifmuth.aZty  bymmztfiic  modeA), 
an  assumption  which  will  be  relaxed  in  Sect.  6. 


Fig.  7.  Radiation  from  an  open-ended  circular  waveguide. 


b’or  a 'I'M  incident,  mode  the  t.urt'aee  current.  incident,  on  the  rim 
no 

a . . 

bar.  only  a '/.—component.*  J_  - «1  *.•,  ; t.ho  d i sport*  ion  ivlat  ion  is: 

J (k  a)  = 0 ; k = BsinO  ; (' 

o t t g 

the  distances  from  Q' , Q"  to  P are  given  by 


r'  = r - asin0cos((J> ' - <f>)  ; r"  - r + asinOcos(<[) ' - <J>)  ; ( 

the  appropriate  elementary  scattering  coefficient  is  (U.3)  with  the  pZui 
sign,  since  surface  currents  have  the  same  sign  at  Q' , Q".  In  conclusion, 
the  total  vector  potential  A = P&  will  be  proportional  to  the  following 
integral: 
d,+7T 

I cos[  3asinGcos((J) ' -<}))] 

I d<D'.  ( 


/-  sin20cos2  (<{> ' - (j))"1-  cos0 


When  Basin©  is  large  the  main  contribution  to  the  integral  (5.3)  comes  from 
integration  points  close  to  <J) ’ = <f>.  When  sinO<<  I/$a  and  3a  is  large,  the 
cosine  term  in  the  denominator  of  (5.3)  does  not  play  a significant  role. 
Accordingly,  we  can  put  (j)'  =(j>  only  in  thd  denominator.  of  (5.3)  and  evaluate 
the  integral  (5-3)  with  the  aid  of  (2.18).  When  all  terms  are  taken  into 
account,  the  final  expression  for  the  vector  potential  is: 

0 

rj  a exp(-iBr)  cos-“  J (BasinO) 

A = - —2— — — exp(U+  - u+)  ( 

-iu>  2r  cos  — cos0  - cosG  ® 

^ 6 

Note  that  the  vector  potential  is  evcAijwh&ie.  finite  in  view  of  the  dispersion 
relation  (5.1). 

We  still  need  to  specify  the  values  of  s and  q which  appear  in  the 
Weinstein  function  U+  (2.7).  It  lias  already  been  noted  in  Sec.  3 that  2trq 


should  equal  the  phase  delay  between  interacting  points  on  the  edge,  such  as 
Q'  and  Q".  In  this  case  of  circular  symmetry,  the  line  .joining  Q'  and  Q" 


intersects  the  axis  of  the  guide,  which  is  a caustic  line.  Accordingly, 


the  usual  phase  delay  of  n/4  should  be  added,  so  that  the  values  of  U—  to 
be  inserted  into  (5.4)  are  the  following: 


+ + , — » 6a  - 11 A 

U (s,  q)  = U (\2&a  cos0,  ) 


(5.5) 


8a  - tt/4 

U (s,  q)  = U (>^28a  cosG  , “ 


:5.6) 


Expression  (5.4)  CO-inCsid&i  with  the  exact  solution  to  the  problem  [20;  , 
with  a minor  modification  in  the  q parameter  of  the  function  U+.  According 
to  Weinstein,  the  proper  value  of  q to  be  inserted  into  (5.5-6)  is: 
q = 4.rg[H^2)(8a)]  (5.7) 


while  in  (5.5-6)  we  have: 

8a  - v 4 
* = — 

Use  of  (2.13)  shows  that  the  difference  between  (5.7)  and  (5.8)  is  of  order 
I/8ir8a. 

The  TE  incident  case  can  be  similarly  treated.  The  surface  current 
no 

on  the  rim  has  only  a (J)  component,  J_  = $;  the  dispersion  relation  is: 

J,(k  a)  = 0 , k = 8sin0  (5.9) 

it  t g 

The  appropriate  elementary  scattering  coefficient  is  (4.2)  with  the  rrUnu i 
sign,  since  surface  currents  have  opposite  signs  at.  Q'  and  Q".  The  elemen- 
tary vector  potential  dA  due  to  the  two  rim  elements  ad<{>  at  Q'  , 0."  can  be 
resolved  into  two  components: 

dA  = dA  cos(4>'  - <{>)$  - dAsin(4>'  -<t>)$  x 9. , (5.10) 

parallel  and  normal  to  respectively.  It  is  readily  seen  that,  the  integral 
of  the  latter  component  is  equal  to  zero,  while  the  integral  of  the  former 


gives  for  t.he  total  vector  potential  A_  = A$  : 


Jin  exp(-ipr)  cos§  Jj(|Ja  sinO) 

A = -it,  exp(ll  - U_)  (‘i.U 

-iu>  2r  cos-^  eosO  - cosO 

2 g 

Note  that  the  vector  potential  is  everywhere  finite  due  to  the  dispersion 
relation  (5.9).  The  function  U should  be  computed  for  the  same  value  (5.8)  of  q. 

Expression  (5.1l)  CO-incidd, 6 with  the  exact  solution  to  the  problem  [20], 
with  the  only  minor  modification: 

9 = - ^ arg[Hj^(Ba)]  (5>12 

Use  of  (2.lU)  shows  that  the  difference  between  (5.12)  and  (5.8)  is  of  order 
3/8irBa. 


CHAPTER  H 


RADIATION  FROM  OPEN-ENDED  CIRCULAR  WAVEGUIDES  CONTINUED: 

AZYMUTHALLY  ASYMMETRIC  MODES 

Radiation  from  an  open-ended  circular  waveguide  is  now  considered  with 

TE  or  TM  excitation,  wherein  the  ^-dependence  of  the  fields  inside  the 
nm  run 

waveguide  is  of  type  exp(im0).  The  analysis  is  more  involved,  although  the 

procedure  parallels  that  for  TE  or  TM  excitation. 

no  no 

First,  TE  mode  excitation  is  considered,  wherein  the  surface  current 
nm 

incident  on  the  rim,  exp(im0),  is: 

mcosQ 

J = J.0+Jz,  6 = - 2 J.  , (6.1) 

— d)  z z „ . .»  <(>  * 

gasin20 

6 

and  the  dispersion  relation  is: 

J'  (k  a)  = 0 , k = BsinQ  . (6.2) 

m t t g 

The  component  of  the  incident  surface  current  produces  a vector  po- 
tential, Aj  = AjZ,  which  can  be  computed  by  using  (U.3),  (5.2)  and  (2.19): 

0 

C.J  a exp(-iBr)  cos—2—  J (BasinG) 

A,  = im  — - exp(im0)  — exp(U+  - U+)  (6.3) 

. G g 

-ito  2r  cos-^r-  cosG  - cosG 

2 6 

where 

. , 3a  - tt/U  - imr/2 

U (s,  q)  = U (^2Ba  cosG,  ) (6.U) 

TT 

Note  the  expression  of  q in  (6. 14 ) . The  term  - tt /1i  accounts  for  the  phase 
delay  at  the  caustic  line  crossing  on  the  axis  of  the  guide.  The  extra  term 
rmr/2  takes  into  account  the  equal  (m  even)  or  opposite  (m  odd)  sign  of  .1 
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at  Q' , Q"  (see  Fig.  7),  as  follows  upon  use  of  (2.8-9). 

A ..A  A 

For  the  vector  potential  = A^<p  - A^tji  x z associated  to  the  component 
J we  have,  upon  use  of  (4.2),  (5-2)  and  (2.20-21): 


i £J,a  exp(-iBr)  cos—  J'(BasinG) 

■im+  tf-  exp(im<{>)  — exp(u‘  - U_)  (6.5) 

gr  ^ 

-iu  2r  cos- ~ — cosG  - cosG 

2 g 


CJ.a  exp(-iBr)  cos^-  m J„( gasinG) 

-im  g—  exp(im<t>)  — exp(U^_  - U ) 

-id)  2r  cos-^-  BasinQ  cosG  - cosG 

C.  6 


(6.6) 


U"(s,  q)  = U 0(iiT  cosG, 


Ba  - rr/4  - mit/2 


(6.7) 


transverse  components  of  the  vector  potential.  Aq  and  A , are  given  by: 


A.  = A’  ; A~  = A'  cosG  - A,  sinQ 

© 2 0 2 / 


(6.8) 


A.  = -i  G 
<P 


.m+J  J a exp(-iBr)  cos| 


J'  (gasinG) 


-id)  2r 


exp(  im<j) ) exp(U  - U ) 

cnsO  - cns0 


cos'-?  cosG  - cosG 

2 g 


J.a  exp(-igr)  2cos— 

. .m © d 

= -i  C -x—  Q- 

-id)  2r  sin_£. 


2c os—  J (gasinG) 

0 exp  ( im<}) ) — A 

sin_£.  sinQ  sinQ 

2 g 


where  U is  given  by  (6.7),  use  has  been  made  of  (6.1)  and  only  terms  up  to 
order  A have  been  retained  in  (6.10). 

Expressions  (6.9-10)  are  now  compared  to  the  expansion  of  exact  results 
[20]  up  to  terms  of  order  A2 included.  The  (^-component  of  the  exact  vector  po- 
tential equals  (6.9)  times 


■J’XX"" 


.’Y 


(m  v 1 + cosO^  cosu  - cosUT 
k^a  J 1 + cosO  2 


osO  - cosO 


provided  that  in  (6.7): 

1 ( , V 

q = - arg  IT  ( 3a ) 
tt  in 

The  O-comoonent  of  the  vector  potential  equals  (6.10)  times: 
/ - A2  = 1 +/-^4cos20 

VV7 


R 


(6.12) 


(6.13) 


Use  of  (2.l6)  shows  that  the  difference  between  the  two  values  of  q is 

small  provided  that  3a  is  even  moderately  large  compared  to  unity.  The 

conclusion  is  drawn  that  uie  ofi  eicmentaiij  Scattering  cee ^ faicient-S  (14.2-3) 

leads  to  good  aeAuits  provided  that  k^a  > m and  the  excitation  is  TE. 

Let  us  now  turn  to  TM  modes,  for  which  the  incident  surface  current 

nm 

J exp(imcj>  > has  only  a z-  component,  J = Jz,and  the  dispersion  equation  is: 

J (k  a)  = 0,  k = 3sin0  (6.1)4) 

m ti  u 

Then,  use  of  (I4.3),  (5.2)  and  (2.19)  will  result  in  a vector  potential 
= Az  given  by  (6.3): 

A = A j (6.15) 

This  vector  potential  will  provide  (Eg,  H^)  fields  only, accord i np  to  (1.6-7), 
while  the  exact  solution  [20]  provides  (Eg,  ) and  ( E^ , Hg)  fields  as  well. 
Comparison  only  between  electric  fields  is  necessary,  since  H = E/t- 

The  exact  expression  for  the  Eg  field,  up  to  terms  A2  included,  is  the 
following : 


.m  + + 

L = -i  7a  2 r exp(im<[>)  exp(U  - U ) 


r 0<r 

cos-^t  J (SasinO) 
2 m 


0 

cos— 


i 11O  - 


m.f  (3asinG)  m 
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cos—  cosO  - cosO 
2 R 


2cos- 


3a 


(6.1b) 


3ns i nO 


(6.17) 


2 o 


where  the  argument  q of  the  function  U is  given  by  : 
q = ~ arg[H^‘‘  ^(3a)l 

As  for  the  TE  case,  difference  in  the  two  values  of  q is  negligible  as  shown 

by  (2.15). 

The  &vu>t  teAm  in  (6.16)  is  exactly  predicted.  The  second  tern,  which 
is  0&  OA.de A (m/Ba)2,  is  not  pAedicted  at  all.  This  term  is  important  only 
for  m = 1 since  , at  variance  with  the  other  one,  is  different  from  zero 
for  0=0.  However,  the  ratio  of  the  first  term  to  the  second  equals: 


for  small  values  of  0,  so  that  the  first  term  dominates  for 


m 0 

0 > tg  — 

23a 

The  exact  expression  for  the  unpAcdicted  ^-component  of  electric 
is,  up  to  terms  of  order  A, 


(6.19) 


field 


exp(-i3r)  <J'  (Basin0 

imCJr^a  exp(im4>)  • m 
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2r 


cos 


(6.20) 


The  conclusion  is  drawn  that,  nOA  TM  excitation,  use  oh  ClemcntaAy  scat- 

nm  ° 

feting  coehhicients  (4.2-3)  leads  to  good  Adults  pAovided  that  3a/m  >>  f; 
discrepancies  with  the  exact  expressions  for  the  fields  arise  only  in  a very 
small  angular  region  centered  at  0 = 0,  for  m = f;  and,  in  the  case  of  linearly 
polarized  excitation,  for  the  cross-polarized  component  of  the  field.  A way 


to  overcome  these  inconveniences 


is  given  in  section  9 • 
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RADIATION  FROM  OTHER  OPEN-ENDED  GUIDING  STRUCTURES 


Some  other  open-ended  guiding  structures  will  now  be  considered  briefly. 
The  vector  potential  associated  to  an  open-ended  coaxial  cable,  of  outer 
and  inner  diameters  2a, 2b  respectively  can  be  easily  computed  using  (4.3) 
and  0^  = 0 for  TEM  excitation.  The  final  result  for  A_  = Az  is  the  following: 


-ico  27rr 


J (3asin0)  - J (BbsinG) 

o o 

sin20 


exp(U  - U ) 
o 


U (s,  q)  = U ( \[3( b-a)Vus0,  — ) 


Note  that  0 is  the  angle  with  the  cable  axis,  I is  the  total  incident  current 
and  only  interaction  between  inner  and  outer  walls  has  been  considered.  When 
3a  <<  J (<6 mall  cable),  eq.  (7.1)  transforms  into: 


A _ £1  exp(-l3r ) g2  s _ M exp ( - i 3r ) , -iMftpD  I , , - 

‘ -iio  8ir2r  4tt  r 2tt 

where  S = ir(b2  - a2  )in  the  aperture  area  and  use  has  been  made  of  (3.18). 

Equation  (7*3)  shows  that  the  cable  aperture  is  equivalent  to  a z-direc- 
ted  electric  dipole  of  moment  ^epsl/2'n . This  is  the  exact  value  which  is  com- 
puted assuming  only  the  electric  field  (twice  as  big  as  the  incident  one)  to 
exist  across  the  aperture  and  then  applying  the  equivalence  theorem  i 3b i . 

For  an  open-ended  (>ViA.p  line  of  width  2b  and  conductor  spacing  a,  the 
E-f  lane  field  can  be  easily  computed  using  (It.-),  (1.6)  and  0 =0  (TEM  exci- 

tation). The  final  result  is  the  following: 


*J| 


21  / . r,  \ 0 sin(0asinO)  , , 

— exp(-ipr)  cos—  ■ . .. expU 


= V 5 


(7.M 


(7 . ^ ) 


U (s,  q)  = U (f 2&a  cosO,  — ) 


(7.6) 


where  I is  the  total  incident  current.  Expressions  (7.^-5)  are  consistent,  in 
the  low  frequency  case,  with  application  of  the  equivalence  theorem,  as  in  the 
case  of  the  coaxial  cable. 

Equations  (7-**-5)  do  not  depend  on  the  strip  width.  Accordingly  they 
are  valid  also  for  an  open-ended  tWO-u)i tAan&miMion  Line. 

The  radiation  from  an  open-ended  rectangular  waveguide  can  also  be  com- 
puted. Results  are  in  good  agreement  with  experimental  data.  They  are  not. 
quoted  hereafter  only  because  some  further  assumption  is  necessary  in  the 
scattering  coefficients  (*4.2-3). 

The  coaxial  cable,  the  strip  and  the  two-wire  transmission  line  are  non- 
dispensive  guiding  structures,  when  excited  under  TEM  conditions.  Accordingly, 
a pulse  is  not  distorted  while  propagating  along  these  guides;  they  are  the 
simplest  suitable  put ie.  fiadidtoU . Transient  radiation  from  these  configu- 
rations has  been  given  in  [26,  36]  neglecting  the  interaction  factor.  Fourier 
inversion  of  this  factor  in  the  time  domain  is  necessary  for  a complete  ana- 


lysis of  transient  radiation  properties  of  the  structures  previously  analyzed. 


3i 


CHAPTER  8 

THE  INCLUSION  OF  A FLARE  ANGLE 


For  the  applications  viewpoint,  it  is  highly  desirable  to  have  scattering 
coefficients  which  can  be  applied  to  fclaAe.d  open-ended  guiding  structures, 
such  as  a horn.  Unfortunately,  no  exact  solution  exists  for  radiation  from 
a flared  parallel-plate  waveguide.  However,  there  is  a way  to  bypass  this 
difficulty. 

Consider  the  sketch  of  Fig.  8 where  the  terminal  section  of  a flared 
parallel-plate  waveguide  is  considered;  let  ct  be  the  flare  angle.  The  field 
distribution  on  the  mouth  z = 0 of  the  flared  guide  will  be  different  from 
that  of  an  unflared  parallel-plate  guide  with  the  same  mouth  dimension  2a 
at  z = 0.  For  small  flare  angles,  however,  this  difference  is  just  a phcu>e. 
eAAOA.  due  to  a change  in  path  length: 


Fig.  8.  Radiation  from  open-ended  flared  parallel-plate  waveguide. 


si  not  cose 


sin2a] 


(1  -h) 


(8.1) 


where  the  only  assuption  used  in  the  last  derivation  is  that  sin2a  <</  . 

Accordingly,  the  phase  error  across  the  mouth  will  be: 

it  y2  & a I sir.  a | 

A4>  = - £2  (1  - - ),  f; 2 = 

2 a2  tt 


(8.2) 


where  B is  the  propagation  constant  in  the  flared  guide  and  f is  a parameter 
describing  the  flare. 

The  field  in  the  half-space  z > 0 can  be  represented  as  the  inverse  Fourier 
transform  (in  wavenumber  space)  of  the  field  distribution  in  the  plane  z = 0. 

We  assume  this  field  distribution  to  be  coincident  with  the  field  distribution 
relative  to  an  unfilaAed  parallel  plate  guide  (of  mouth  dimension  2a)  tvnoA  a 


pkaAe  eM.cn: 


A<J)  = ^ £ 2 U - 2)  , -a  < x < a 

A<J>  = 0 | x I > a 


(8.3) 


(8.M 


Accordingly,  the  field  for  z > 0 and,  in  particular,  the  far  field,  will  be 
obtained  as  the  convolution  of  the  (known)  field  radiated  by  the  (unflared) 
parallel  plate  waveguide  and  the  F ounleA  tnanA  faonm  the  complex  exponential 
06  (8.3-1*) : 

-a  a 

ip  = J exp(  iuy  )dy  + ^ J exp(iuy)  exp[i  77  ( 1 - ) ] dy 


exp( iuy )dy 


where  u is  the  y-eomponenl.  oi'  spatial  finquoiicv  :m.i  oqua  , s ps- i in|>  (jiv  Fig.  .) 
ill  i ■ i •: i I space. 

Computation  of  (8.S)  is  rather  involved  and  requires  t.he  use  of 
The  final  result  is  the  following: 

= S(u)  - ^ j — - jTj-  |f(^-  - C)  exp(iua)  - f(~  + K ) exp(-iua)J 

- j r |g(^|-  -£)  exp(iua)  - g(~  +f.)  exp  (-iua)]|  (8.6) 


When  C is  small  eq.  (8.6)  can  he  simplified  by  expanding  the  auxiliary 
functions  f and  g in  the  neighbors  of  ua/TTt) . Using  (2.26)  we  get: 


ip  = 6(u) 


PL  £. 

IT 


i ua  „,ua,  . ua  , uaH 

aJ  L T Ts*j 


Mu!  . vU)  (3.7) 


As  a check,  it  can  be  noted  that,  for  £-*  0,  the  last  bracketed  term  in  (8.7) 
reduces  to  -i(TTh/ua)2,  upon  use  of  (2.28).  Accordingly,  the  second  term  of 
(8.7)  is  proportional  to  £2  and  ip  reduces  to  a 6 - function  for  £ = 0,  as  it 
should  be. 

The  new  vector  potential  A,  relative  to  the  flared  parallel  plate  wave- 


guide, is  given  by: 


■ n 


' ) A(u  - u’ )du 1 = A ( u ) + 


TOO 

/ 


)A(u  - u' )du' 


( • 


where  the  vector  potential  A relative  to  the  unflared  parallel  plate  system  is 
given  by  (3.1-2)  and  (3.8-9). 

Due  to  the  discontinuous  behavi  >r  ; . 4-2  ) f the  fund  i n f(l  and  • • . 
• 1-  range  f the  Integral  which  appears  should  he  split  into  the  two 

parts  (-»,  0)  and  (0,  <*>).  For  the  latter  we  have,  upon  repeated  integration 
by  parts  and  use  of  (2.26): 


Now,  the  bracketed  term  approaches  -i  and  remain?  constant,  a?  t > uf,  as  follows 
upon  use  of  (2.28);  there  is  no  serious  error  in  taking  this  term  approximately 
constant.  Then,  integration  by  parts  of  the  first  integral  at  right  hand  side 
of  (8.9)  gives: 
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dA 


<Mu'  )A(u  - u')du'  = +i£‘  -f~  + iL 

0 du 


oo  _ 

ip2  f L - — t— 

' J [A  t dt  tdt  t " 


d A sint 
dt  t 


dt  (8.10) 


Now  the  rapidly  varying  part  of  A(ua  - t)  is  of  type  exp(+it),  when  the  aperture 
is  large.  Accordingly,  d2A/dt2  ~ -A  and  (8.10)  transforms  into: 


f. 


( u 1 )A(u  - u'  )du ' = +iC2  ^ 

du 


+ l?f 

*/0 


;in  t - t cos  t , , , 

—3 A(ua  - t.  )dt 


(8.11) 


When  also  the  remaining  range  of  the  integral  (8.8)  is  taken  into  account,  we 
get: 


OO 

L ^(u')A(u  - u ' )du ' = + 2if,2 
■'U  du 


,.r  2 dA 
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3in  t - t cos  t 

7~1 


|^A  ( UE 


ua  - t. ) 


- A(  ua  + t) 


dt 


(3.12) 


Eq.  (8.12)  provides  a way  for  computing  the  scattering  coefficients  analogous 
to  (h.2-3)  and  appropriate  t.o  the  open  ended  flareu  waveguide. 
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CONCLUDING  KKMARKG 


A new  formulation  has  been  presented  for  radiation  from  open-ended 
guiding  structures.  It  is  essentially  a field  theory  which  calculates  the 
field  as  an  integral  over  the  rim  of  the  guide.  No  problems  of  ray  tracing, 
or  boundaries,  or  caustics  arise.  The  theory  can  even  be  extended  to  open 
guides  with  a terminal  flare  angle.  Although  presented  only  for  the  front 
space,  the  formulation  can  be  extended  to  the  rear  space. 

The  key  point  of  the  theory  is  the  canonical  problem  of  a parallel -plate 
waveguide.  It  is  not  unexpected  that  when  the  theory  is  applied  to  structures 
far  from  the  canonical  one,  some  discrepancies  are  found.  This  is,  for  example, 
the  case  for  a circular  waveguide  excited  by  azymuthally  asymmetric  modes, 
where  the  incident  field  on  the  rim  of  the  guide  presents  a ^-variation  with 
no  counterpart  in  the  original  canonical  problem  of  a parallel  plate  waveguide. 
Obviously,  we  can  take  as  baiic  canonical  pAobCem  that  ofj  the  cinculan  waveguide 
with  azijmuthatt  ij  aSijmme.tAdc  mode  i;  and  derive  mete  general,  although  more  compli- 
cated, catteiing  coefifiicienti  which  would  replace  those  of  eqs.  (a.  1;. 

such  a way  the  applicability  of  this  new  formulati  n woul  i certainly  be 
enlarged  at  the  expense  of  some  formal  eompl i cations . 

Guiding  structures  which  are  worth  analysing  are  the  corrugated  circular 
waveguide  and  the  rectangular  waveguide. 

For  the  former,  the  different  type  of  boundary  conditions  requires  the 

use  of  incident,  eiectfiic  as  well  as  magnetic  surface  current:-.  A first  rough 

% 

way  to  handle  the  problem  would  be  the  use  of  scattering  coefficients  for  the 
magnetic  currents  obtained  from  those  of  the  electric  one  by  duality.  A 
better  approach  would  be  to  use  the  solution  of  scattering  by  a half-plane  with 


two  different  surface  impedances  (in  the  particularly  suitable  form  [37])  to 
model  new  scattering  coefficients. 

For  the  latter,  the  formulation  just  presented  for  the  open-ended  parallel 
plate  waveguide  is  already  able  to  give  good  results.  However,  the  use  of 
di^eA&nt  dispersion  relations  for  the  E and  H planes  is  necessary.  This  is 
certainly  a rather  unsatisfactory  point  and  is  again  due  to  the  necessity  of 
handling  a structure  too  different  from  the  canonical  one. 

We  believe  that  the  limitations  of  this  theory  are  not  intrinsic  but 
rather  reflect  its  early  stages.  It  is  probably  worth  developing  further  this 
novel  formulation,  which  presents  a simplicity  and  a number  of  attractive  fea- 
tures not  available  in  current  theories  of  high  frequency  radiation. 
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